Influence of spin waves on transport through a quantum-dot spin valve 
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We study the influence of spin waves on transport through a single-level quantum dot weakly 
coupled to ferromagnetic electrodes with noncollinear magnetizations. Side peaks appear in the 
differential conductance due to emission and absorption of spin waves. We, furthermore, investigate 
the nonequilibrium magnon distributions generated in the source and drain lead. In addition, we 
show how magnon-assisted tunneling can generate a fully spin-polarized current without an applied 
transport voltage. We discuss the influence of spin waves on the current noise. Finally, we show 
how the magnonic contributions to the exchange field can be detected in the finite-frequency Fano 
factor. 
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I. INTRODUCTION 

Quantum dots coupled to ferromagnetic electrodes 
have generated much interest in the recent past due to 
their possible applications in spintronics. Experimen- 
tally realizations of such systems include small metal- 
lic grains, 1-3 semiconductor quantum dots, single-wall 
carbon nanotubes, 10-14 as well as single molecules 15 con- 
tacted by ferromagnetic electrodes. From a theoretical 
point of view, quantum-dot spin valves, i.e., quantum 
dots coupled to noncollinearly magnetized leads are of 
particular interest for a number of reasons. First, they 
show a nonequilibrium spin accumulation on the quan- 
tum dot due to spin-dependent tunneling. This spin ac- 
cumulation has the tendency to block transport through 
the system and therefore gives rise to a tunnel magne- 
toresistance. 16-20 On the other hand, there is a preces- 
sion of the dot spin due to an effective exchange field 
acting on the dot. 21-24 This exchange field arises due to 
spin-dependent virtual tunneling processes between the 
dot and the leads. It helps lifting the spin blockade by 
allowing the dot spin to precess out of its blocking po- 
sition. It is, thus, responsible for a deviation from the 
harmonic behavior of the conductance as a function of 
the angle between the magnetization directions of the 
leads. 21 In the nonlinear transport regime, the interplay 
between spin accumulation and spin precession leads to 
a broad region of negative differential conductance. 22 A 
more direct access to the exchange field is provided by 
the finite-frequency Fano factor that shows a resonance 
signal at the Larmor frequency associated with the ex- 
change field. 24 A splitting of the Kondo resonance due 
to the exchange field was predicted using the numerical 
renormalization group 25-29 and observed in recent exper- 
iments. 13,15 ' 30 Recently, a new way to access the exchange 
field in a quantum-dot spin valve with an additional su- 
perconducting electrode was proposed. 31 

Further theoretical studies addressed the current 
noise 32-35 and the full-counting statistics 36 of single-level 
quantum dots coupled to ferromagnetic electrodes, as 



well as transport through multi-level dots, 35 ' 3 double 
dots, 20,38 ' 39 and carbon nanotube dots 34,40-43 coupled to 
ferromagnetic electrodes. 

While the transport properties of quantum-dot spin 
valves have now been investigated in quite some detail, 
more extended models that include, e.g., the possibility 
to excite spin waves in the ferromagnetic electrodes have 
not yet been addressed. In this work, we consider the 
influence of spin waves on transport through a quantum- 
dot spin valve using a suitable extension of the real-time 
diagrammatic transport theory. 44-47 On the one hand we 
want to analyze the deviations from the idealized sys- 
tem. To this end, we study the modifications of the con- 
ductance which we find to be particularly pronounced 
for large polarizations of the leads. Here, the magnonic 
side peaks can exhibit negative differential conductance 
but can also surmount the ordinary conductance peaks 
depending on the magnetic configuration. We, further- 
more, show that the excitation of spin waves can lead 
to an increased as well as to a decreased Fano factor. 
Additionally, we demonstrate how the magnonic modifi- 
cations to the exchange field can be detected in the finite- 
frequency noise. On the other hand, we want to address 
the question if the spin waves can generate completely 
new effects. To this end, we analyze the nonequilibrium 
distribution of the magnons which we find to be different 
for the source and drain electrode. Furthermore, we show 
how the magnons can drive a completely spin-polarized 
current without any external bias voltage. 

The paper is organized as follows. In Sec. II we in- 
troduce our model. The extension of the real-time dia- 
grammatic technique to systems containing spin waves is 
presented in Sec. III. Our results are presented in Sec. IV. 
We conclude by giving a summary in Sec. V. 



II. MODEL 

We consider transport through a single-level quantum 
dot weakly coupled to ferromagnetic leads with non- 
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FIG. 1. (Color online) Sketch of a quantum-dot spin valve 
with spin wave excitations in the leads and coordinate system 
used. 



collinear magnetizations via tunneling barriers. In or- 
der to describe spin waves which may be excited in the 
leads, we model the lead magnetizations as macroscopi- 
cally large spins localized in the respective lead as shown 
in Fig. 1. 

The Hamiltonian of our system is the sum of five parts 



H = H, 



dot 



H r + Ht, 



coupl • 



(2.1) 



The first term, H^ot = J2a £c\c a + U c^c^^c^ describes 
the quantum dot in terms of a single, spin-degenerate 
level with energy e measured relative to the Fermi ener- 
gies of the leads in equilibrium and Coulomb energy U 
for double occupation. 

The ferromagnetic leads are modeled as reservoirs of 
itinerant electrons, H r = J2rka=± £ r^a\^ (y a r \ ia , r = L,R. 
Here, a rka denotes the creation operator for electrons in 
lead r with momentum k and spin a which have energy 
e r k- The spin quantization axis of each lead is chosen 
parallel to the direction of its magnetization e r . 

Due to the noncollinear geometry, it is convenient to 
quantize the spin on the dot in the direction perpendicu- 
lar to the magnetizations of the leads, cf. the coordinate 
frame indicated in Fig. 1. In this case, the tunnel Hamil- 
tonian describing the coupling between dot and leads is 
given by 



-fftun — 



rk 



V2 



•Vk-I 



+4- (-e 40 " /2 ct + e-^/ 2 c ; )] + h.c, (2.2) 

where <f) r denotes the angle enclosed between the mag- 
netization of lead r and the x axis. It should be noted 
that though individual terms violate spin conservation, 
the total tunnel Hamiltonian is spin conserving. 

Using the Holstein-Primakoff representation, the local- 



ized spins are expressed in terms of bosonic operators 

S rz =S-btb r , (2.3) 
S r+ = (J2S^bJb~r \ b r , (2.4) 

(2.5) 



S r - = bt f \/2S~btb r 



Since the magnetizations are macroscopic quantities, 
their fluctuations will be small and hence b\b r <C S. 
We, therefore, can restrict ourselves to the leading-order 
terms when expanding Eqs. (2.3)-(2.5) in 1/5*. The spins 
are then treated as a bosonic degree of freedom with en- 
ergy u) b , H spin = ^2 r bj b b\b ri corresponding to optical 
magnons with a single, dispersionless mode. The cou- 
pling between the localized spins and the spins of the 
itinerant electrons in the leads is given by 



H, 



coupl 



— J ^ ■ Sf> 

r 

- - P ( a rk+ a rk+ - <4k- a rk- ) (2.6) 

rk 

+ V2S (aj k+ a rk -& r + a,^ k _a rk+ 6 r 



The terms in the second line are independent of the 
magnon number. They can be absorbed into the part 
of the Hamiltonian describing the lead electrons, yield- 
ing spin-dependent electron energies, £ r ko- = e-rk — crJS. 
As a consequence, the densities of states will also ac- 
quire a spin dependence, p r cr (w). In the following, we 
will refer to the electrons with larger (smaller) density 
of states at the Fermi energy as majority, a = + (mi- 
nority, a = — ) spin electrons. Since in transport only 
a small energy window at the Fermi energy is relevant, 
we will take the densities of state to be independent of 
energy, p r<rT — p r<(T (Ep). In this case, the tunnel ma- 
trix elements t r can be related to the tunneling rate r r ± 
for a majority /minority spin electron from lead r via 



tunneling rate r r = 



We furthermore define the average 

(r r+ + r. r _)/2. 

The asymmetry between majority and minority spins 
is then characterized by the spin polarization p r = (p,. i+ — 
Pr,—)/(Pr,+ + Pr.-)- We note that in general for a non- 
constant density of states, the total number of majority 
spin electrons may be larger or smaller than the total 
number of minority spin electrons, depending on the pre- 
cise form of the band structure, the size of the splitting 
and the position of the Fermi energy. Hence, the local- 
ized spin may point parallel or antiparallel to the major- 
ity spin direction corresponding to positive or negative 
values of the polarization p r . We note that a noncon- 
stant density of states gives rise to a modified exchange 
field, Eq. (3.8), see below. 

The terms in the last line of Eq. (2.6) describe spin-flip 
interactions between lead electrons and the magnons. In 
order to remove this interaction, we apply to the Hamilto- 
nian the canonical transformation 48 H = e A He~ A with 
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generator 

A = - ^ A ( a rk+ a rk-^ - <4k- a r-k+&r) (2-7) 
rk 

where A = Jv / 2S/(wf,+e r k+-£rk-) = J^2S /(w h -2JS). 
Neglecting terms of order A 2 that can be absorbed into 
the energies e r ker, the transformed Hamiltonian takes the 
form 

H = ffdot + H r + H spin + H tun (2-8) 
with the transformed tunneling Hamiltonian 

Htun = ^2 ~h ( fl rk+Crt + a rk-Cri) + h.C. (2.9) 
rk 

where 

<H = (l - Y bf r b r) (f r '\ + ^ l2 C^ (2.10) 

Sri, = (l - Y b r b r) (-^'\ + e^C^ (2.11) 

+ Xb r + e-^c^j . 

In writing down the transformed tunnel Hamiltonian, 
Eq. (2.9), we neglected terms of order A 2 that do not 
involve bosonic operators as these terms do not yield a 
contribution to order A 2 in the diagrammatic expansion. 
The canonical transformation gives rise to new processes 
in the tunnel Hamiltonian. Apart from tunneling events 
that are already present for the ordinary quantum-dot 
spin valve, 22 we encounter additional terms in which a 
magnon is emitted/absorbed and the spin of the tunnel- 
ing electron is flipped. Hence, the total tunnel Hamilto- 
nian does not conserve the electron spin any longer but 
just the sum of electron spin and angular momentum of 
the magnons. 



III. REAL-TIME DIAGRAMMATIC 
TECHNIQUE 

In order to investigate the transport properties of our 
system, we extend the real-time diagrammatic technique 
developed in Refs. 44-47 and adapted to systems with 
ferromagnetic leads in Refs. 21 and 22 to include the 
spin-wave degrees of freedom. 



A. Reduced density matrix and master equation 

The basic idea of the real-time diagrammatic technique 
is to integrate out the non-interacting, fermionic degrees 



of freedom of the leads in order to obtain an effective de- 
scription of the reduced system which is characterized by 
the state of the quantum dot and the number of magnons 
in the left and right lead. The quantum dot can be either 
empty, occupied with a spin up or a spin down electron, 
or doubly occupied. We denote these states as |0), | f), 
| J,), and \d) with energies E — 0, E^ = E^ = e, and 
Ed = 2s + U, respectively. The number of magnons is 
characterized by |n) = \ni,,ns}. The total energy of a 
state |£) = |x, n) is then given by Eg = E x + (nL,+nR)u)b- 
In order to keep the dimension of the Hilbert space finite, 
we introduce a maximal magnon number iV max for actual 
computations and check that our results are independent 
of the cut-off value. For the parameters chosen in the 
analysis below, it turned out that it is sufficient to take 
into account at most four magnons in each lead. 

After tracing out the fermionic degrees of freedom in 
the leads, the system is described by a reduced density 
matrix p rcd with matrix elements = (£2|p red |£i)- For 
the diagonal matrix elements we introduce the abbrevia- 
tion P e = P|. 

In the stationary state, the density matrix elements 
obey a generalized master equation of the form 

o = f£ = l (s i2 - Eil ) P ^ + j2 Ki p i- 

Here, the kernels wf 2 *? describe transitions due to the 
dot-lead coupling. Similarly as for the diagonal den- 
sity matrix elements, we introduce the short-hand form 
W||, = W^'. The kernels that enter the master equa- 
tion are defined as irreducible self-energy diagrams for 
the dot propagator on the Keldysh contour and can be 
expanded perturbativcly in the tunnel coupling strength 
r. In the following, we will restrict ourselves to terms 
that are of first order in T. The diagrammatic rules nec- 
essary for the evaluation of the diagrams are summarized 
in Appendix A. 

In the following, we will always assume the magnon 
energies to be much larger than the tunnel coupling, 
LJb 3> r. In this case, the master equation for the matrix 
elements which are off-diagonal in the magnon number 
becomes i(E^ 2 — £Jj 1 )P| i 2 = to first order in the tunnel 
coupling. This implies that only matrix elements diago- 
nal in the boson number have to be taken into account. 
In this case, the reduced density matrix takes a block- 
diagonal form p rcd = pdot <£> Pboson, where each block can 
be written as 

/4 \ 

o Pu Pi o 

Pf n P U 

V o o o p dn J 

To allow for a compact notation, we write the den- 
sity matrix as a vector consisting of blocks of the form 

{Pon, Pw P -i!n' ^tn^ The normalization of the 

reduced density matrix can then be cast into the form 



e T p red = 1, where e T is a vector consisting of blocks of 
the form (1,1,1,1,0,0). 

To allow an easier physical interpretation of the density 
matrix elements, it is convenient to express them in terms 
of the average dot occupations Pon, -Pin = Pfn + Pin an d 
Pdn and the average spin on the quantum dot 



where 



dS n 
~dt 



e e e 4nm^< mer ' 



r m x 



dS n 

~dT 



(3.5) 



(3.6) 



P 



R 



tn 



p\ _ pi p p 

rt _ - 4.il fn Q _ fn ^4.n 



Z Z i. 

(3.2) 

in the presence of n magnons. The set of master equa- 
tions can then be split into one describing the average 
charge and one describing the average spin on the dot. 

The master equation for the occupation probabilities 
is given by 



dt 



^M (r) P, +V {r) S 

X' 



(3.3) 

(r) 

Here, , m denotes transition rates from state \ m to 

(r) 

state xn, while V^m characterizes the dependence of the 
dot occupation on the accumulated spin. Their precise 
form is given in Appendix B. The sum over m only gives 
contributions if in a process involving the left (right) lead 
the number of magnons in the right (left) lead is kept 
fixed while the number of magnons in the left (right) 
lead changes by at most one due to the conservation of 
angular momentum in each tunneling event. 

This restriction in the number of excited magnons is 
different from the case of a quantum dot coupled to a 
vibrational degree of freedom. In the latter case, the 
number of phonons that can be excited in a tunneling 
event is limited only by the applied bias voltage while 
the number of phonons that can be absorbed is only 
limited by the number of excited phonons, 49 giving rise 
to a large number of conductance sidebands 50-52 as well 
as to a phonon distribution width that is nonperturba- 
tive in the electron-phonon coupling. 49 Furthermore, the 
transition rates in the vibrational case show a nontriv- 
ial dependence on the initial and final number through 
the Franck-Condon factors that influence the transport 
properties crucially, e.g., by leading to a suppression of 
transport for small bias voltages. 51 ' 52 

The dot spin obeys the Bloch-type equation 



dS n 

dt 



dSn 

dt 



dSn 

dt 



rcl 



dS n 

dt 



(3.4) 
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FIG. 2. (Color online) Exchange field as a function of the 
level position e for p r = 0.9 (upper panel) and p r = 0.1 (lower 
panel). Other parameters are U = 50kBT, ujb = WksT. 



dS n 

dt 



= s b x£b« 



(3.7) 



prcc 



The dynamics of the dot spin is governed by three terms. 
The first one, Eq. (3.5), which depends only on the occu- 
pation probabilities, describes the accumulation of spin 
on the dot due to electrons tunneling onto the empty dot 
or electrons leaving the doubly occupied dot. Here, the 
sum over m is subject to the same restriction as in the 
master equation for the occupation probabilities. The 
second term, Eq. (3.6), which is proportional to the ac- 
cumulated spin, describes the decay of the dot spin due 
to tunneling out of electrons or tunneling in with a spin 
opposite to the dot spin, forming a spin singlet on the 
dot. The precise form of the functions F^ m and 
that enter the accumulation and relaxation term is given 
in Appendix B. Finally, the third term, Eq. (3.7), de- 
scribes the precession in the exchange field generated by 
virtual tunneling between the dot and the leads. It is 
given by 



prcc 
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e r — < (1 - X 2 n r ) Pr [$ r (e) - $ r (e + 17)] + tt 

7T Z 



21n^ (1 + p r )(l + n r )$ r (e - u h ) 



+(1 -p r )n r $ r (e + io b ) + (1 +p r )n r $ r (e + U - u b ) - (1 -p r )(l + n r )<& r (£ + JJ + o^) 



(3.8) 



where $ r (a;) 



Re 



-Mr) 



and ^ is the 



2 1 " 2tt 

digamma function. Compared to the ordinary quantum- 
dot spin valve, 22 the exchange field contains new terms 
proportional to A 2 which arise form virtual tunneling 
processes that emit or absorb magnons in the interme- 
diate state. Furthermore, the terms already present in 
the absence of spin waves experience a magnon-number- 
dependent renormalization 1 — X 2 n r that has the ten- 
dency to reduce the strength of the exchange field. In 
Fig. 2 we show the exchange field as a function of the 
level position. The new terms give rise to side peaks 
and dips. As the original exchange field is rather small 
for small polarizations, the magnonic features tend to be 
more pronounced for small polarizations. 

The logarithmic divergency in Eq. (3.8) that is cut- 
off by the bandwidth of the lead electrons W (in the 
following, we assume W — 100/cbT) arises as the rate for 
emitting and absorbing a magnon differ from each other. 
If the system is, e.g., in a state without magnons and 
the dot is singly occupied, only a spin down electron can 
leave the dot to the leads by a magnon- assisted process 
while only a spin up electron can enter the dot in such a 
process. Hence, only the energy of the spin-down state is 
renormalized by the magnonic processes, thereby giving 
rise to a diverging energy shift between spin up and spin 
down electrons and hence to a diverging exchange field. 

In order to take into account the finite life time of spin 
waves, e.g., due to scattering from phonons or electrons 
or due to magnon-magnon interactions, we include phe- 
nomenological relaxation terms — - (P xn — Pn q ICn -^x") 
and — - (S x „ — P° q J2 a S X n) into the master Eqs. (3.3) 
and (3.4). These terms describes a relaxation to- 
wards the equilibrium distribution of magnons P® q = 

e -^b/k B T B _ e -u: b /k B T B ^ Qn a time gcale T _ We al j ow 

for the most general case of a magnon temperature Tb 
that differs from the electron temperature T. 

Besides making our model more realistic, the relax- 
ation terms also ensure that the magnon number in the 
drain lead remains finite. While the emission of a magnon 
when tunneling out of the dot transfers the electron from 
minority to majority spin and thus gives rise to a rate 
proportional to 1 + p r , the absorption flips the spin in 
the opposite direction and therefore yields a rate propor- 
tional to 1 — p r . Hence, without any relaxation mecha- 
nism, the number of magnons would grow without bound 
in the drain lead. 



B. Current and current noise 

We define the current through the tunnel barrier r as 
the change in the number of electrons in lead r multiplied 

with the electron charge, I r = — ie[H, Ylka a rkcr a rkcr]- 
While in the stationary state the current at zero fre- 
quency obeys (I) = (7l) = —{Ir) due to current conser- 
vation, at finite frequencies the current measured in the 
source-drain circuit in general differs from the currents 
through the tunnel barriers due to occurrence of displace- 
ment currents. These can be taken into account accord- 
ing to the Ramo-Shockley theorem 53-55 by considering 
the total current as the sum of the current I r through 
each tunnel barrier weighted with the corresponding ca- 
pacitance C r of the barrier, I — (Cl/l + Cr/r)/(Cl + 
Cr). As the capacitances of the tunnel barriers are much 
less sensitive to the geometry than the tunnel couplings, 
in the following, we assume symmetric capacitances while 
allowing for asymmetric tunnel couplings. We therefore 
have I = (7l + 7r)/2. Diagrammatically, the current can 
be evaluated as 



{I) 



2h 



-e 1 W J P, 



(3.9) 



where P denotes the vector of density matrix elements 
and are kernels in which one internal vertex is re- 
placed by a current vertex. Since the current operator 
equals the tunnel Hamiltonian apart from constant fac- 
tors, these replacements only give rise to factors of ±1 
depending on the position of the vertex on the contour 
and whether the electron tunnels in or out of the dot (cf. 
Appendix A for details). 

We define the frequency-dependent current noise 2 as 
the Fourier transform of S(t) = {I(t)I(Q)) + {I(Q)I(t)) - 

S(w) = dt ((/(t)7(0)) + (J(0)7(t))) (e^ + e~^) 

-4tt6(u){I). (3.10) 

By choosing a symmetrized expression for the current 
noise, we restrict ourselves to real noise measurable in 
a classical detector 56 in contrast to the unsymmetrized 
expression which is complex and describes detector- 
dependent emission and absorption processes. 57-62 

Following Ref. 24, we evaluate the finite-frequency 
noise diagrammatically to first order in the tunnel cou- 
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FIG. 3. (Color online) Transport processes that arise in a 
quantum-dot spin valve with spin wave degrees of freedom in 
the leads as well as their dependence on the polarization p. 



pling as 



S(w) = ^e T [W J/ + W 1 (Uq 1 ^) - W) W 1 ] P 

- 2ir5(uj)(i) 2 + (lo -> -u). (3.11) 



Here, the matrices W 7/ correspond to self-energies with 
two tunnel vertices replaced by current vertices. Further- 
more, we introduced the free propagator of the quantum 
dot, 



™Xix' 1 "X2X' 2 



£ X2 - £ X i - w + * 0+ 



IV. RESULTS 



(3.12) 



In this section, we discuss our results for the current, 
conductance and current noise of the quantum-dot spin 
valve in the presence of spin waves. Unless stated oth- 
erwise, we assume a symmetric system with equal po- 
larizations for both leads, p\, = PR = P, and equal 
tunnel couplings, Tl = Tr = T/2. Furthermore, we 
assume that the bias voltage is applied symmetrically, 
V h = -Vk = V/2, too. 
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FIG. 4. (Color online) Schematic of the differential conduc- 
tance in the e — V plane. The labels at the conductance lines 
indicate whether a magnon is emitted (E) or absorbed (A) in 
the source (S) or drain (D) electrode. The dotted horizontal 
line marks the position of the cuts shown in Fig. 5. 



quantum dot. Hence, the rate for this process is given by 
A 2 ?t,l(1 +p)T/4. Similarly, a minority spin electron can 
emit a magnon to become a majority spin that tunnels 
onto the dot. Here, the corresponding rate is given by 
A 2 (l + n L )(l-p)r/4. 

Furthermore, there are two processes involving the 
magnons in the drain lead. Now, a minority spin of the 
drain can leave to the dot, emit a magnon and end up 
as a majority spin. The rate for this process is given by 
A 2 (l + n R )(l +p)T/4. The opposite process which starts 
with a majority spin, absorbs a magnon and ends up as a 
minority spin finally has a rate given by A 2 tir(1 — p)T/A. 

Hence, we find that for a given spin direction on the dot 
the magnonic processes have the opposite dependence on 
the polarization as the normal processes. We will show 
in the following how this unconventional polarization de- 
pendence gives rise to a number of interesting transport 
properties. 



B. Magnon-assisted tunneling 



A. Transport processes 

Before discussing the transport properties in detail, we 
give an overview of the different transport processes that 
can occur in the presence of spin waves in the leads. In 
Fig. 3, we show all transport processes that can arise. 
As for the ordinary quantum-dot spin valve, majority 
(minority) spin electrons can tunnel from the source onto 
the dot as well as from the dot into the drain lead. The 
corresponding rates are given by (1 ±p)F/4. 

Additionally, we now have processes that involve the 
emission/absorption of magnons. For example, a major- 
ity spin from the source can flip its spin by absorbing a 
magnon to become a minority spin that ends up on the 



We now turn to the discussion of the differential con- 
ductance G = dl/dV for arbitrary bias and gate volt- 
ages. In Fig. 4, we schematically show the conductance 
for collinear magnetizations 63 . In addition to the reso- 
nances of the dot level with the Fermi energy of the left 
and right lead which are marked by the thick black lines, 
a number of sidebands related to the emission and ab- 
sorption of magnons occur. At the red (light) dashed 
lines labeled ES, the tunneling electrons become able 
to emit magnons in the source lead, while at the blue 
(dark) dashed lines ED the emission of magnons in the 
drain lead becomes energetically possible. The dashed 
green (dark) lines indicate where transport through the 
dot becomes possible by absorbing magnons to enter the 
empty dot with both level above the two Fermi energies 
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(AS) or to leave the doubly occupied dot with both lev- 
els below the two Fermi energies (AD). Furthermore, 
there are additional peaks marked by the orange (light) 
dotted lines AS and AD at which it becomes possible to 
absorb a magnon in order to enter the doubly occupied 
state which otherwise would have been out of reach en- 
ergetically. Conductance sidebands inside the Coulomb 
blockade regime are absent because there a spin accumu- 
lates on the dot in such a way as to suppress all magnonic 
processes. 

As discussed above, the different magnonic processes 
depend on the polarization of the leads as 1+p or 1— p de- 
pending on whether they involve a majority or minority 
spin electron in the lead. As a consequence, the differ- 
ent sidebands will have different strengths. In particular, 
while line ES is suppressed with increasing the polariza- 
tion, line ED is increased. Hence, the conductance map 
no longer exhibits particle-hole symmetry — £ due 

to the presence of the magnons. 

As discussed in Sec. II the polarization of the leads can 
take either sign: a positive (negative) sign indicates that 
the majority spins of the carriers at the Fermi energy are 
parallel (antiparallel) to the direction of the macroscopic 
magnetic moment. We now discuss how the conductance 
plot changes for the different combinations of positive 
and negative polarizations in the two leads. The case of 
two negative polarizations is related to the case of two 
positive ones by the transformation e -+ U — s. This 
maps all lines involving the emission of magnons onto 
lines involving the absorption of magnons and vice versa 
and therefore exchanges the factors 1+p and 1 — p in 
the magnonic rates. For polarizations of opposite signs, 
the conductance plot becomes particle-hole symmetric. 
However, now the strength of the magnonic side bands 
depends on the direction of current flow. While they are 
suppressed for one direction by 1 — p they are enhanced 
by 1 + p in the other direction. In the rest of this work, 
we will restrict ourselves to the case in which both po- 
larizations are positive. All other cases can be related to 
this one by the above symmetry considerations. 

We now turn to the discussion of the conductance at 
fixed bias voltage shown in Fig. 5, corresponding to a cut 
along the dotted horizontal line in Fig. 4. In the upper 
panel of Fig. 5, we show the differential conductance for 
parallel magnetizations as a function of the level position 
for a given bias voltage. For small polarizations, there are 
four large peaks at which the dot levels are in resonance 
with either the left or right Fermi energy, corresponding 
to the thick black lines in Fig. 4. These are accompanied 
by much smaller conductance peaks at distances ±cjf,. 
From left to right, these correspond to lines AD, ED, 
ES, ED, ES, and AS in Fig. 4. The lines correspond- 
ing to the absorption of magnons (AS and AD) are only 
present for sufficiently high magnon temperatures, oth- 
erwise the number of magnons in the leads is too small 
to make these side peaks visible. 

For large polarizations, there are two new effects. On 
the one hand, the two side peaks associated with lines 
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FIG. 5. (Color online) Differential conductance as a function 
of the level position for V = 30fc B T, U = 50fc B T, ujb = 
10fc B T, T L = r R , r = 2/r, A = 0.3, T B = 5T for small 
and large polarizations. In the upper panel, we show the case 
of parallel magnetizations, while in the lower panel they are 
chosen to be antiparallel. As in Fig. 4, the labels of the side 
peaks indicate the emission (E) and absorption(E) of magnons 
in source (S) and drain (D) lead. 



labeled ES are suppressed by 1 — p and become invisi- 
ble. On the other hand, we find that the side peak corre- 
sponding to left line ED in Fig. 4 shows a strong negative 
differential conductance. This is due to the formation of 
a trapping state. A spin down electron can leave the dou- 
bly occupied dot by exciting a magnon in the drain lead 
leaving the dot in the spin up state. This blocks further 
transport as tunneling in of an electron is suppressed by 
the small density of states for spin down electrons in the 
source lead. 

We now turn to the case of antiparallel magnetizations. 
For small polarizations, the conductance resembles the 
one for the parallel configuration because the effects of 
the finite polarization are only weak. For large polar- 
izations, there is a large spin accumulation on the dot 
which is antiparallel to the magnetization of the drain 
lead. When transport takes place through the states |0) 
and |cr) only, the rate for normal tunneling into the drain 
is suppressed by 1 — p while the rate for tunneling via a 
spin flip is proportional to A 2 (l +p). Hence, for large po- 
larizations and not too small couplings to the magnons, 
the side peaks can dominate over the main peaks. On 
the other hand, for transport through the states \a) and 
\d), the bottleneck is given by the tunneling in of elec- 
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FIG. 6. (Color online) Difference between the average number 
of magnons in the source and drain lead in the parallel (upper 
panel) and antiparallel (lower panel) configuration for r = 
100/r, p = 0.9 and Tb = T. Other parameters are the same 
as in Fig. 5. 



tron. This bottleneck cannot be overcome by magnonic 
processes. Hence, we do not find large side peaks in this 
case. 

In summary, this means that for large polarizations 
the magnons tend to decrease the current for transport 
through singly and doubly occupied dot in the parallel 
configuration while they tend to increase the current for 
transport through the empty and singly occupied dot for 
antiparallel magnetizations. As a consequence, the tun- 
nel magnetoresistance will be reduced by the magnons 
for all transport regimes. 



C. Nonequilibrium magnon distribution 

As the tunneling processes in the quantum-dot spin 
valve can emit and absorb magnons, they will give rise 
to a nonequilibrium magnon distribution in the source 
and drain lead. As we discussed in Sec. IV A, the rate 
for absorbing a magnon in the source lead is enhanced 
compared to the rate for emitting a magnon. Hence, 
the average magnon number in the source will be re- 
duced compared to the equilibrium distribution. For the 
drain lead, the situation is reversed. Here, the rate for 
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FIG. 7. (Color online) Current at zero bias voltage vs. level 
position for U = 50fc B T, u = 20fc B T, T L = T K , t = 1/T, X = 
1/3, Tb = 10T at different polarizations. Inset: Increased 
relaxation time leads to a reduction of the magnon-driven 
current (shown here the peak value at e — WksT) since it 
allows a stronger cooling of the magnons. 



emitting a magnon is enhanced resulting in an increased 
magnon number compared to equilibrium. In Fig. 6 we 
plot the difference between the average magnon number 
in the source and drain lead as a function of gate and bias 
voltage. We find that indeed the number of magnons in 
the drain lead is larger than the one in the source lead. 
For parallel magnetizations, this effect occurs for all gate 
voltages, while for antiparallel configurations the effect 
occurs only for transport through the states |0) and \a) 
as only in this case the magnons can help to overcome 
the spin blockade on the dot. As they do this very effi- 
ciently, the nonequilibrium effects are more pronounced 
here compared to the parallel configuration giving rise to 
a larger deviation between the average magnon number 
in source and drain. 



D. Magnon-driven electron transport 

In the following, we show how an asymmetric cou- 
pling to the magnons in the left and right electrode can 
lead to a completely spin-polarized, magnon-driven cur- 
rent at zero bias voltage. A finite current quite gener- 
ically requires a breaking of detailed balance, W^> — 
e p(E t -E t ,)yy^^ whii e usually detailed balance is bro- 
ken by a finite applied bias voltage, here we break de- 
tailed balance by different equilibrium temperatures for 
the magnons and electrons, Tb ^ T. Similar effects occur 
quite generally when the system couples asymmetrically 
to external fields and breaks detailed balance. 46,64 Exper- 
imentally, this mechanism has been realized by coupling 
microwaves to the gate electrode defining a quantum dot 
in a two-dimensional electron gas 65-68 and in a carbon 
nanotube 69 giving rise to photon-assisted tunneling. 

For the quantum-dot spin valve, there are different 
ways to achieve an asymmetric coupling to the magnonic 
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degrees of freedom. We can either choose asymmetric 
tunnel couplings, Tl 7^ Tr, different couplings to the 
spin waves X r or different polarizations p r . In the fol- 
lowing, we choose a system with one ferromagnetic and 
one normal lead, i.e., we have magnons only in the left 
lead and pn = 0. In order to violate detailed balance, 
we choose the equilibrium distribution of the magnons 
to have a larger temperature than the electrons in the 
leads. We note that to describe magnon-assisted tunnel- 
ing, we do not have to keep track of the magnon numbers 
in the leads explicitly. Instead, we could also integrate 
out the magnons and described them as an additional 
bath with temperature Tb. 

In Fig. 7, we show the magnon-assisted current for van- 
ishing transport voltage as a function of level position 
for different polarizations. For < e < ujb, a completely 
spin-polarized current flows from the ferromagnet into 
the normal lead. A spin down electron in the ferromag- 
net can absorb a magnon and tunnel onto the empty 
dot. Subsequently, the electron can leave the dot either 
into the ferromagnet or into the normal metal giving rise 
to a net current into the normal metal. Increasing the 
polarization enhances the probability of the spin-flip pro- 
cesses that populate the quantum dot. Furthermore, the 
tunneling rates for the spin-flipped electron back into 
the ferromagnet are decreased. Hence, with increasing 
the polarization of the ferromagnet we find an increased 
magnon-assisted current. 

Similarly, for — U — ujb < e < —U, a completely spin- 
polarized current flows from the normal lead into the 
ferromagnet. In this case, a spin up electron can leave 
the quantum dot into the ferromagnet by absorbing a 
magnon. Afterward, electrons from the left as well as 
the right lead can tunnel onto the dot, yielding a net 
current from the normal to the ferromagnetic lead. Since 
the rate for the absorption processes now is proportional 
to 1 — p, these processes become strongly suppressed for 
large polarizations. 

In addition to these current plateaus whose width 
scales with the magnon energy, there are additional peaks 
at e = and e = — U whose width is given by the elec- 
tron temperature. They occur since at resonance with 
the leads, the system may flip the spin on the quantum 
dot which tends to align as to inhibit magnon-absorption 
processes, thereby making the processes described above 
possible again and yielding a finite current. 

As just illuminated, the key ingredient to the zero-bias 
current at Tb > T is the absorption of magnons (simi- 
larly, for Tb < T it would be the emission of magnons). 
This explains why the height of the current plateaus 
scales linearly with the average number of magnons. Fur- 
thermore, it explains why the current is reduced when the 
relaxation time is increased (see inset of Fig. 7). In this 
case, the system absorbs magnons in the processes dis- 
cussed above. Since it takes more time to relax to the 
equilibrium magnon distribution now, the average num- 
ber of magnons is reduced and therefore also the current 
is reduced. 



If we compare our results for the magnon-driven cur- 
rent to the case of a current driven by photon-assisted 
tunneling, we find that although the basic mechanism of 
absorbing bosons to gain energy is the same in both cases, 
there are nevertheless some striking differences. First, 
in contrast to the photon case, there are no processes 
which involve the absorption of more than one magnon 
which is a result of angular momentum conservation dur- 
ing tunneling, cf. the discussion in Sec. III. Second, while 
the photons drive an unpolarized current, the magnons 
yield a fully spin-polarized current since in absorption 
processes they couple only to minority spins. This also 
leads to a breaking of particle-hole symmetry for finite 
values of pl- 



E. Current noise 

So far we discussed the current and conductance of a 
quantum-dot spin valve with spin wave degrees of free- 
dom. We now turn to the discussion of the current 
noise which can provide additional information about the 
transport processes occurring in the system. 



1. Zero- frequency noise 

We start with the discussion of the zero-frequency 
noise. As shown in Fig. 8, where we plot the zero- 
frequency Fano factor F(w = 0) = S(oj = 0)/(2eJ) 
as a function of the level position, for sufficiently high 
polarizations, a quantum dot coupled to ferromagnetic 
leads with parallel magnetizations will exhibit super- 
Poissonian Fano factors in the regime where only two 
charge states of the dot contribute to transport. This is 
due to a dynamical spin blockade, 24,32 ' 33,70 ~ 75 where mi- 
nority spin electrons block transport due to their longer 
dwell time on the dot. Thereby, they chop the current 
into bunches of majority spin electrons. 

If the excitation of spin waves is energetically possible, 
the Fano factor still remains super-Poissonian. However, 
it is significantly reduced when transport takes place 
through the states |0) and \a) (right plateau in Fig. 8) 
while it is slightly increased for transport through the 
states \a) and \d) (left plateau in Fig. 8). The reduction 
of the Fano factor on the right plateau can be understood 
qualitatively by taking into account that the emission of a 
magnon in the drain lead gives the possibility for a minor- 
ity spin to leave the dot. Therefore, magnonic processes 
reduce the waiting time between bunches of majority spin 
electrons and reduce the Fano factor in consequence. On 
the left plateau, due to the magnonic processes, there are 
now two ways to get into the blocking state. Either a spin 
down electron tunnels to the drain in a normal tunneling 
event or it does so in a spin-flip process by exciting a 
spin wave. Hence, the bunching effect becomes stronger 
in this region and the Fano factor is increased. 
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In order to gain a more quantitative understanding of 
the influence of magnons on the Fano factor, we define 
an effective polarization of the leads which takes into 
account the magnonic processes. In the absence of spin 
waves, the ratio between the rate for transferring a spin- 
up electron from the dot into the lead and the rate for 
transferring any electron from the dot to the lead equals 
(1 +p)/2. In the presence of spin waves, there are new 
processes (cf. the discussion in Sec. IV A) which have a 
different dependence on the polarizations. We therefore 
define effective polarizations that take into account the 
presence of the new magnonic processes via 



FIG. 8. (Color online) Fano factor as a function of the level 
position for parallel (upper panel) and antiparallel (lower 
panel) magnetizations. In both cases the Fano factor can be 
increased or decreased by the coupling to the magnons. Pa- 
rameters are V = 40fc B T, U = 50k B T, uj b = Wk B T, T L = T R , 
p = 0.9, t = 2/r and T B = T. 



1 + PL.cff _ W^f"L"R,0"L"R + ^t"L + lnR,Qn L n R 

+ W tnL + lnR,0riLnR ~t~ ^4<nL — l^R,0nLriR 

1 + PR,cff _ Wp nijnR ^ n - LnR + Wora L raR-l,traL"R 

2 IForiLnn^nLnR + ^0n L n B 4n L nR + "07JL7JR-l,tnL«R + WonL«R + l,i"LnR 



(4.1) 
(4.2) 



r 



In the above expressions, we set the magnon numbers 
that occur in the rates equal to the average magnon num- 
ber found from the solution of the master equation. This 
is a reasonable approximation for the parameters chosen 
here as the system has a very high probability to be in a 
state with zero magnons. 

Hence, the system with spin waves can be interpreted 
as a quantum dot coupled to ferromagnetic leads with dif- 
ferent polarizations pl and pr (Note that here different 
signs of pl and p R correspond to antiparallel magneti- 
zations and are not related to the negative polarizations 
discussed in Sec. II). For such a system, the Fano factor 
can be computed analytically and is given by 



F(0) = 



(1 - 2p LPR + p| )( 5 + 2 PLPR + pl) 



(3 - 2p L p R 



Pl) 



(4.3) 



on the right plateau. On the left plateau, the same ex- 



pression holds with pl and pji exchanged. In Fig. 9, 
we show the Fano factor as a function of pa for various 
values of pi, in both transport regions. It is interesting 
to note that for large polarizations, the Fano factor de- 
pends quite sensitively on the precise polarization value 
such that the spin waves can significantly alter the Fano 
factor. Furthermore, the behavior on the two plateaus is 
completely different. 

Using the picture of effective polarizations, we find that 
for the parameters of Fig. 8 where due to Wb ^ Tb and 
t ~ 1/r the average number of magnons is nearly van- 
ishing, we have PL,cff = 0.90 and p RyC s — 0.83. According 
to Eq. (4.3), this yields F(0) = 2.1 for the right plateau 
in good agreement with the numerically obtained value. 
Similarly, by exchanging pl and p R in Eq. (4.3), we find 
F(0) = 4.8 for the left plateau, which is also in good 
agreement with the full calculation. 
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FIG. 9. (Color online) Zero-frequency Fano factor for a quan- 
tum dot symmetrically coupled to collinear magnetized leads 
as a function of the polarization of the right lead for different 
values of the polarization of the left lead. In the upper panel, 
transport takes place through the empty and singly occupied 
dot only, while in the lower panel, it takes place through the 
singly and doubly occupied dot only. 



In the antiparallel case, shown in the lower panel of 
Fig. 8, the zero-frequency Fano factor tends to unity for 
large polarizations without magnons. In this case, an 
electron leaving the dot to the drain lead will be easily 
replaced by another majority electron from the source 
lead which then blocks transport for a long time. Hence, 
the tunneling out events become uncorrelated and the 
Fano factor Poissonian. As in the case of parallel mag- 
netizations, the presence of magnons will result in an 
effective polarization of the drain smaller than the true 
polarization. As can be read off from Fig. 9, this results 
in a reduced, sub-Poissonian Fano factor for the right 
plateau, while it leads to an increased, slightly super- 
Poissonian Fano factor for the left plateau. For trans- 
port at the right plateau, the magnons allow the spin up 
electron on the quantum dot to leave to the drain lead 
by flipping its spin. This reduces the waiting times be- 
tween tunneling events to the drain and therefore makes 
the Fano factor sub-Poissonian. On the other hand, at 
the left plateau, after a spin-up electron has left the dot 
by emitting a magnon in the drain, another spin up elec- 
tron will enter the dot. In the next step, this scheme will 
either be repeated or a spin down electron will leave the 
dot, resulting in the spin-blockaded state. Hence, the 
magnons can initiate minibunches that lead to slightly 
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FIG. 10. (Color online) Frequency-dependent Fano factor for 
perpendicular magnetizations for V — 25fc B r (upper panel) 
and V — 75fc B T (lower panel). Other parameters are U = 
50fc B T, e = 10fc B T, oj B = 10fc B T, T L = 2r R , r = l/r L , 
p = 0.8, T B = T. 

super-Poissonian Fano factors. 

Since for perpendicular magnetizations, where the 
Coulomb plateaus become modulated due to the influ- 
ence of the exchange field, 24 one does not find any new 
effects beside the already discussed reduction of the Fano 
factor for the right and an increase for the left plateau, 
we do not discuss this case here further. 



2. Finite-frequency noise 

We finally turn to the discussion of the finite- frequency 
Fano factor F(uj) for frequencies comparable to the tun- 
nel couplings, uj < r. As was discussed in Ref. 24, in 
contrast to the average current which is only sensitive to 
the average spin on the dot, the finite-frequency current 
noise provides more direct access to the spin dynamics on 
the quantum dot. In the following, we discuss how the 
frequency-dependent Fano factor can be used to detect 
the magnonic contributions to the exchange field. The 
exchange field can only give rise to a precession of an 
accumulated spin if the exchange field and spin are not 
collinear to each other. Since the exchange-field contribu- 
tion of the each lead is parallel to its magnetization direc- 
tion, the leads magnetization should not be collinear. To 
be specific, we choose perpendicularly magnetized leads 
in the following. 
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In Fig. 10, we show the frequency-dependent Fano fac- 
tor for two different bias voltages and different coupling 
strengths to the spin waves. For a vanishing coupling to 
the magnons, the Fano factor exhibits a peak at the Lar- 
mor frequency associated with the exchange field. As 
the strength of the coupling to the spin waves is in- 
creased, the height of the resonance peak is reduced. Fur- 
thermore, the position of the resonance peak is shifted 
as the strength of the exchange field is altered. Since 
the magnonic contributions to the exchange field have 
a different bias dependence as the standard contribu- 
tions, 22 the magnitude and sign of the Larmor frequency 
shift depends on the bias voltage. Hence, measuring the 
frequency-dependent Fano factor as a function of bias 
voltage can provide experimental access to the magnonic 
exchange field contributions. 

V. CONCLUSIONS 

We investigated the influence of spin waves excited in 
the electrodes on transport through a quantum-dot spin 
valve. We found that the excitation and discexcitation 
of spin waves gives rise to conductance sidebands whose 
strength depends on the polarization and the relative ori- 
entation of the magnetizations. We, furthermore, found 
that the transport through the system gives rise to a 
nonequilibrium occupation of the magnons with an in- 
creased magnon number in the drain and a decreased 
magnon number in the source lead. For a system that 
couples asymmetrically to the magnons in the source and 
drain, we showed how magnon- assisted tunneling leads to 
a completely spin-polarized current without an applied 
bias voltage. Finally, we studied the current noise. We 
found that the magnons can increase as well as decrease 
the Fano factor depending on the dot states involved in 
transport. Additionally, we showed that the frequency- 
dependent Fano factor provides access to the magnonic 
contributions to the exchange field. 
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Appendix A: Diagrammatic rules 

The diagrammatic rules for computing the kernels 
W^ 2 |, 2 in frequency space are given by 



1. Draw all topological different diagrams with tun- 
neling lines connecting vertices on either the same 
or opposite propagators. Assign to the four cor- 
ners and all propagators states |x, ^l, n R) and cor- 
responding energies E x +LOb(n-^ + 7t.r) as well as an 
energy u for every tunneling line. 

2. For each time interval on the real axis confined by 
two adjacent vertices, assign a resolvent 1/(AE + 
irj) where AE is the difference between left and 
right-going tunneling lines and propagators. 

3. For each tunneling line involving lead r, the di- 
agram acquires a factor where the sign 
on the Fermi function depends on whether the line 
runs forward (— ) or backward (+) with respect to 
te Keldysh contour. 

4. For each pair of vertices connected by a 
tunneling line the diagram is multiplied by 

^<£|2rtl£a)<3l4tl&) + V<CMU(ai^l6> 

where |£ a ) and \£' a ) (|£&) and \£' b }) are the states 
that enter and leave the vertex the tunneling line 
begins (ends) at, respectively. The operators cfij 
are defined in Eq. (2.10) and (2.11). In evaluating 
the above matrix elements, take into account only 
terms up to order A 2 . 

5. Assign a factor of (— l) a + h where a is the num- 
ber of vertices on the lower propagator and b is the 
number of crossings of tunneling lines. 

6. Sum over all leads r. 

7. Integrate over all energies of tunneling lines. In 
the sequential-tunneling regime only one tunneling 
line is involved. In this case, the frequency inte- 
gration reduces to a simple application of Cauchy's 
formula. 

8. For the computation of W J and W 11 , one, respec- 
tively two tunnel vertices are replaced by current 
vertices. These give rise to a factor +1/2 if they are 
on the upper (lower) branch of the Keldysh contour 
and describe the tunneling of an electron into the 
right (left) or out off the left (right) lead. Otherwise 
they result in a factor —1/2. 



Appendix B: Master equation 

In this appendix, we give expressions for the various 
functions that enter the master equations (3.3) and (3.4). 
Introducing the abbreviations 
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r r = (l+p r )n r f+(e -oj b ) + (1 - p r )(l + n r )/+(e + 
s. r = (1 + p r )(l + n r )f~{e + U — Ub)-(1- p r )n r f~(e + w b ), 
x r = f-(e) + f+(e + U), 
y r = f-(e)-f+(e + U), 
z r± = (1 +p r )(l + u r )/~(e - w&) ± (1 -p r )n r f~(s + w 6 ) 

±(1 + Pr )n r f+{e + U-ut,) + {l- p r )(l + n r )f+(e + U + w 6 ), 

we can write the quantities M^nx'm as matrices in the basis |0), |1), |c?): 

2(l-A 2 n,,)/+(e)-A 2 r r (1 - A 2 n r )/" (e) 
M«=r r | 2(l-A 2 n r )/+(e) -(1 - A 2 n r )y r + ^z r+ 2(1 - A 2 n r )/" (e + 17) 

(1-A 2 n r )/+ (e + U) -2(1- \ 2 n r )f-(e + U)- X 2 s r 

V ' / rVf/-(e-w 6 ) 

' rc r I 2T r _f+(e + uj b ) 2r r+ /"( £ + U - u b ) 

r r _/+(£ + C7 + w 6 ) 



M 



(r) 



(r) 



2 

A 2 



r r _/-(e + w 6 ) 

(l + n r ) \ 2T r+ f+{e-u Jb ) 2r r _/-(e + f/ + u b ) 

r r+ /+( £ + C/-o; 6 ) 



(Bl) 
(B2) 
(B3) 



I 

(r) (r) 

The vectors V^n and i^ nm that enter the master equa- and 
tion for the occupations, Eq. (3.3), and spin, Eq. (3.4), 
respectively, can be written as vectors in the basis |0), 
|1) and \d) as 



(r) _ 



2p r T r (l — \ 2 n r ) —y r + \ 2 z r - 
V fr{e + U) 

r r _/~(£ + uj b ) 
o 

T r+ f-(e + U -u b ) 

T r+ f-(e-u> b ) 

VA-i = -AV I 

T r _f-(E + U + uj b ) 



Vi%=\ 2 (l + n 




F<;)= Pr r r (i-AV)| _2_y L p^_ 

T T+ f+{e-u b ) 



i 7 . 



(■'■) 



nn+l 



?{r) 



A 2 

y 

A 2 



(l + nr) I 

T r _f-(e + U + uj b ) 

r r _f+(e + u b ) 




r r+ /-(e + ^-w 6 ) 



Finally, we have 



Q( r ) = (i — X 2 n r )T r x r + —Y r z r+ . 



(B7) 
,(B8) 
(B9) 
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